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Abstract. We consider the partial transpose of the spin reduced density matrix of two 
disjoint blocks in spin chains admitting a representation in terms of free fermions, such as 
XY chains. We exploit the solution of the model in terms of Majorana fermions and show 
that such partial transpose in the spin variables is a linear combination of four Gaussian 
fermionic operators. This representation allows to explicitly construct and evaluate the 
integer moments of the partial transpose. We numerically study critical XX and Ising 
chains and we show that the asymptotic results for large blocks agree with conformal 
field theory predictions if corrections to the scaling are properly taken into account. 
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1. Introduction 

In the last decade, the understanding of the entanglement content of extended quantum 
systems boosted an intense research activity at the boundary between condensed matter, 
quantum information and quantum held theory (see e.g. Refs, pj as reviews). The 
bipartite entanglement for an extended system in a pure state is measured by the well- 
known entanglement entropy which is the von Neumann entropy corresponding to the 
reduced density matrix of one of the two parts. One of the most remarkable results in 
this held is the logarithmic divergence of the entanglement with the subsystem size in 
the case when the low energy properties of the extended critical quantum systems are 
described by a 1+1 dimensional conformal invariant theory |21 El HIE!. 

Conversely, when an extended quantum system is in a mixed state (or one considers 
a tripartition of a pure state and is interested in the relative entanglement between two 
of the three parts) the quantification of the entanglement is much more complicated. A 
very useful concept is that of partial transposition. Indeed it has been shown that the 
presence of entanglement in a bipartite mixed state is related to occurrence of negative 
eigenvalues in the spectrum of the partial transpose of the density matrix [6j. This 
led to the proposal of the negativity [7] (or the logarithmic negativity) which was later 
shown to be an entanglement monotone [8j, i.e. a good entanglement measure from a 
quantum information perspective. Compared to other entanglement measurements for 
mixed states, the negativity has the important property of being easily calculable for an 
arbitrary quantum state once its density matrix is known (and indeed for this reason it 
has been named a “computable measure of entanglement” 0 )- 

Recently, a systematic path integral approach to construct the partial transpose 
of the reduced density matrix has been developed and from this the negativity in 1+1 
dimensional relativistic quantum held theories is obtained via a replica trick [9] - This 
approach has been successfully applied to the study of one-dimensional conformal held 
theories (CFT) in the ground state [9111 (Jj. in thermal state [TTJ|T23, and in non-equilibrium 
protocols mmmm, as well as to topological systems [16, Ej. The CFT predictions 
have been tested for several models uni nn uni uni uni iso] , especially against exact results 
EDI HU H31 EH for free bosonic systems (such as the harmonic chain). Indeed for free 
bosonic models, the partial transposition corresponds to a time-reversal operation leading 
to a partially transposed reduced density matrix which is Gaussian [[221 and that can be 
straightforwardly diagonalised by correlation matrix techniques |23| 2H [25]. It should 
be also mentioned that there exist some earlier results for the negativity in many body 
systems 

[2HE21EHE31EZIESIE31EQIEI]- 

In the case of free fermionic systems (such as the tight-binding model and XY 
spin chains) the calculation of the negativity is instead much more involved. Indeed 
the partial transpose of the reduced density matrix is not a Gaussian operator and 
standard techniques based on the correlation matrix [25] cannot be applied. In view of the 
importance that exact calculations for free fermionic systems played in the understanding 
of the entanglement entropy [3j [32, 33] [M, 331 [36], it is highly desirable to have an exact 
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representation of the negativity also for free fermionic systems. A major step in this 
direction has been very recently achieved by Eisler and Zimboras [3TJ who showed that 
the partial transpose is a linear combination of two Gaussian operators. Unfortunately, 
it is still not possible to extract the spectrum of the partial transpose and hence the 
negativity, but at least one can access to integer powers of the partial transpose which 
are the main ingredient for the replica approach to negativity. 

In Ref. [37] only truly fermionic systems have been considered and not spin chains 
that can be mapped to a fermionic system by means of a (non-local) Jordan-Wigner 
transformation. Indeed, in the very interesting case of two disjoint blocks in a spin chain 
the density matrix of spins and fermions are not equal [39] sni SI] and this consequently 
affects also the partial transposition, as already pointed out in [57], In this manuscript 
we fill this gap by giving an exact representation of the partial transpose of the reduced 
density matrix for two disjoint blocks in the XY spin chain and from this we calculate 
the traces of its integer powers. These turn out to converge to the CFT predictions in 
the limit of large intervals. 

The manuscript is organised as follows. In Sec. [2] we describe the model and the 
definition of the quantities we will study. In Sec. [3] we review the results of Ref. [01] for 
the moments of the spin reduced density matrix of two disjoint blocks. In Sec. [4] we move 
to the core of this manuscript deriving an explicit representation of the partial transpose 
of the spin reduced density matrix as a sum of four Gaussian fermionic matrices. This 
allows to obtain explicit representations for the moments of the partial transpose. In Sec. 
[5] we use the above results to numerically calculate these moments up to n — 5 for the 
critical Ising model and XX chain and carefully compare them with CFT predictions by 
taking into account corrections to the scaling. I 11 Sec. [ 6 ] we numerically evaluate and 
study the moments of the partial transpose for two disjoint blocks of fermions and again 
we compare with new CFT predictions. Finally in Sec. [7] we draw our conclusions. I 11 
appendix [A] we report all the CFT results which we needed in this manuscript. 


2. The model and the quantities of interest 


In this manuscript we consider the XY spin chains with Hamiltonian 

L 


3 = 1 


Hxy = - 5 Y ( + h 


3 j+i 


( 1 ) 


where erf are the Pauli matrices at the j-th site and we assume periodic boundary 
conditions cr £ +1 = cr“. For 7 = 1 Eq. (Jl| reduces to the Hamiltonian of the Ising model 
in a transverse field while for 7 = 0 to the one of the XX spin chain. The Hamiltonian (|T|) 
is a paradigmatic model for quantum phase transitions [38]. In fact, it depends on two 
parameters: the transverse magnetic field h and the anisotropy parameter 7. The system 
is critical for h = 1 and any 7 with a transition that belongs to the Ising universality 
class. It is also critical for 7 = 0 and \h\ <1 with a continuum limit given by a free 
compactihed boson. 
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Figure 1. We consider the entanglement between two disjoint spin blocks A\ and A 2 
embedded in a spin chain of arbitrary length. The reminder of the system is denoted by 
B which is also composed of two disconnected pieces Bi and B 2 . 


The Jordan-Wigner transformation 



( 2 ) 


maps the spin variables into anti-commuting fermionic ones {cj,cj} = 
these fermionic variables the Hamiltonian (|TJ) becomes 


Sij. In terms of 



(3) 


where we neglected boundary and additive terms. This Hamiltonian is quadratic in the 
fermionic operators and hence can be straightforwardly diagonalised in momentum space 
by means of a Bogoliubov transformation. 

For the study of the reduced density matrices it is very useful to introduce the 
Majorana fermions |B] 



(4) 


which satisfy the anti-commutation relations {a r ,a s } = 25 rs . 

2.1. Quantities of interest 

The main goal of this manuscript is to determine the entanglement between two disjoint 
intervals in the XY spin chain. We consider the geometry depicted in Fig. [l] a spin chain 
is divided in two parts A and B and each of them is composed of disconnected pieces. 
We denote by Ai and A 2 the two blocks in A = Ai U A 2 , Bi is the block in B separating 
them, while B 2 is the remainder. 

The reduced density matrix of A is p A = Tr#p = Tr£|T)(T|, where we are mainly 
interested in the case in which |T) is the ground state of the XY chain, even if the 
results of this paper apply to more general cases such as excited states, non-equilibrium 
configurations, finite temperature etc. The bipartite entanglement between A and B is 
given by the well-known entanglement entropy 


(5) 


S A = -Tr p A In p A , 

or equivalently by the Renyi entropies 



( 6 ) 


which in the limit n —> 1 reduce to the entanglement entropy, but provide more 
information since it is related to the full spectrum of p A m- 
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In the case of two disjoint blocks in the XY spin chain, the Renyi entropies for integer 
n (or equivalently the moments of the reduced density matrix p/f) have been explicitly 
constructed in Ref. m- However it is still not possible to find the analytic continuation 
to arbitrary complex values of n and consequently the entanglement entropy. This is very 
similar to the CFT counterpart where also one can calculate only integer moments of the 
reduced density matrices (see Appendix [A] for a summary of the CFT results of interest 
for this paper). 

However, we are here interested in the entanglement between A\ and A 2 . A measure 
of this entanglement is provided by the logarithmic negativity defined as follows. Let us 
denote by e,- 1 ^) and \ef^) two arbitrary bases in the Hilbert spaces corresponding to A\ 
and A 2 . The partial transpose of pa with respect to A 2 degrees of freedom is defined as 


(1) (2), T 2 , (1) (2)v _ /„(1)„(2) | n m i„(i)„(2)\ 


(e- (j |p 


e i) = 


'\pM 


(7) 


and then the logarithmic negativity as 

£ = In Up? 11 =lnTr|pJ|, (8) 

where the trace norm \p T f \ | is the sum of the absolute values of the eigenvalues of p T f. 

A systematic method to compute the negativity in quantum held theories has been 
developed in Ref. PHD] and it is again based on a replica trick from the integer moments 
of the partial transpose 

Tr (p5T, (9) 


which turn out to admit different analytic continuations from even and odd n (usually 
denoted as n e and n Q respectively). Consequently, the logarithmic negativity is given by 
the replica limit 


£ = lim InTr (pj)" 6 , 

n e —s>l 


( 10 ) 


performed on the even sequence of moments. Unfortunately, in the case of two disjoint 
intervals, it is very difficult to perform this analytic continuation (see Appendix 0- 
although the integer moments Tr(/rJ ) n are known for the most relevant conformal held 
theories. It is however possible to extract very useful information about entanglement 
and about the partial transpose already from the knowledge of the moments PHD]. 

Finally, in Refs. PHD] the ratios 


R 


n 


Tr(p?) 

Tr Pa 


( 11 ) 


have been introduced because of some 


cancellations (see Appendix [A]) and since anyhow 


£ = lim ln(R n J, 

n e —>l 

given that Trp^ = 1. 


( 12 ) 
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In this section we review the results of Ref. mi. where the Renyi entropies of two disjoint 
blocks A = Ai U A 2 for the XY spin chains have been computed. These results are the 
main ingredients to construct the integer powers of the partial transpose which will be 
derived in the following section. We will denote the number of spins in A 1 and A 2 with 
l\ and l 2 respectively and the remainder of the system B contains a region separating Ai 
and A 2 denoted as B i, as pictorially depicted in Fig. [TJ 

In a general spin 1/2 chain the reduced density matrix p a = Tr^|T) (T| of A = A 1 UA 2 
can be computed by summing all the operators in A as follows [3] 

^ = ^h'E(U a ijU a ?’ < 13 ) 

Vj j&A jeA 


where j is the index labelling the lattice sites and Uj € { 0 , 1 , 2 ,3}, with a 0 = 1 the identity 
cr 2 = u v and a 3 = a z the Pauli matrices. The multipoints correlators 


matrix and a 1 = a x , a 2 = cr y and cr 3 = a 


in (13) are very difficult to compute, unless there is a representation of the state in terms 
of free fermions. 

For the single interval case, the Jordan-Wigner string Y\ m< j a m Eq. (| 2 j) maps the 
first i spins into the first t fermions [3j so the spin and fermionic entropy are the same. 
In general, calculating a fermionic density matrix is very easy. Indeed by Wick theorem 
they assume a Gaussian form: 

ex P (1 E r . s a rW rs a s ) 


Pw = 


Tr[exp Q a r W rs a s )] ’ 


(14) 


with W a complex antisymmetric matrix. This density matrix is univocally identified by 
the correlation matrix 


T rs Tr (Qj r pwa s ') d rs , 


(15) 


which (in matrix form) satisfies T = tanh(lF/2). 

Unfortunately, the same reasoning does not apply to two disjoint blocks because 
the fermions in the interval B\ separating the two blocks contribute to the spin reduced 
density matrix of Ai U A 2 [3Sj SD, E]- In particular, the following string of Majorana 
operators appears in a crucial wajj|] 

Pb i = JJ (ia 2 j -ia 2j ). (16) 

j'eBi 


Similarly, we also introduce the strings of Majorana operators Pa, and Pa 2 , defined as in 
Eq. (16) but taking the product over j e A\ and j G A 2 respectively. These operators 
satisfies Pfj 1 = Pc, for any C = A 1 , A 2j B x . 


Moving back to the computation of the spin reduced density matrix (13) in terms of 
fermions, we first notice that, since the XY Hamiltonian commutes with ]/[. a z - (where j 
runs over the whole chain), only the expectation values of operators containing an even 


t To compare with Ref. [4T], set (Si) there = {PaA here and (*S f )there — (-PBi) here* 
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number of fermions are non vanishing. Thus, the numbers of fermions are either even or 
odd in both Ai and A 2 . This leads us to decompose the spin reduced density matrix 
of A = Ai U A 2 as [CL] 

PA Peven T Pb i Podd) (17) 

with 

^{ 0 , 0 2 > 0 * 0 ;, °l°\’ < 18 > 


Peven — 


2h+h 1, Podd - 2h+f , 2 

even odd 


where we introduced the ‘shorts’ 0\ and 0 2 for arbitrary products of Majorana fermions 
in Ai and A 2 respectively, and the notation ^j) even Eodd) means that the sum is restricted 
to operators Oi and 0 2 containing an even (odd) number of fermions. It is also convenient 
to rewrite as 

_ 1 + Pb i 1 — Pr x 

Pa — -x-P+ + 


P-, 


P± = Pe 


i Podd; 


(19) 


2 2 

where 1 is the identity matrix and (1 ± PgJ/2 are orthogonal projectors. Moreover, 
the matrices p± are unitary equivalent (indeed Pa 2 P±Pa 2 — P- f) and commute with Pb 1 
because they do not contain Majorana fermions in B\. As a consequence, we have that 


Pa = 


1 + 


Pl + 


B i 


P-. 


( 20 ) 


Taking the trace of (20) and employing that Pa 2 P±Pa 2 = p-f, one hnds 

Tr p\ = Tr p" ± . (21) 

The matrices p± are fermionic but they are not Gaussian, i.e. they are not proportional 
to the exponential of a quadratic form, as we will discuss below. 

At this point we are ready to consider the ground state of the XY chain with density 
matrix pw = |T)(T|, whose correlation matrix T is given by Eq. (15). The fermionic 
reduced density matrix of A = A± U A 2 is 

1 £< 0l 0 2 }0*0i, (22) 


Pa = 


2 ^i 


even 

odd 


where we recall that (0\0 2 ) = 0 when the numbers of fermionic operators in 0\ and 0 2 
have different parity. In order to take into account the effect of the string P Bl defined in 
Eq. (16), it is useful to introduce the auxiliary density matrix 


B, _ Tr fl (P Bl |W)(*|) 
Pa — 


{Pb,) 


in which the normalisation Trp^ 1 = 1 holds. By using that 


Pa 2 Pa 2 


J e A 2 , 
3 ^ ^ 2 ) 


one hnds 


Pmp\Pm = 


2 ^i 

1 

2 ^i+& 


even 

odd 


Y,(Oi0 2 )OlO\ - ^ ^(0,0 2 ) 0*0j, 


(23) 

(24) 

(25) 

(26) 


odd 


even 
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where p 2 is the number of Majorana operators occurring in 0 2 . Then, from Eqs. (22) 
and (26) it is straightforward to get that the density matrices (18) become 

Pa 2 Pa l Pa 2 


Peven 


Pa + Pa 2 p\Pa 2 


Podd — (P Bl ) 


Pa 1 


(27) 


2 ' ' 2 
Plugging (27) into ( [T7| ), one finds that the spin reduced density matrix is a linear 
combination of four fermionic Gaussian operators. Since these operators do not commute, 
they cannot be diagonalised simultaneously and therefore we cannot find the eigenvalues of 
the spin reduced density matrix that would give the entanglement entropy. Nevertheless, 
Tr p n A for integer n can be computed through Eq. (21) by providing the product rules 
between the four Gaussian operators occurring in (27) in terms of the corresponding 
correlation matrices, that are denoted by 


ri — 


r 2 = r 


Pa 2 Pa p a 2 ' 


r 3 = r 


Pa 


1 5 


r “ “ I R,."' r V 


(28) 


where is the correlation matrix of a Gaussian density matrix p as in Eq. (15). 
Obviously, Yi is the fermionic correlation matrix, i.e. the one of the free fermions without 
the Jordan-Wigner string (studied in detail in Ref. [43i|). 

Following Ref. |4T], we can introduce the restricted correlation matrix to two 
fermionic sets/blocks C and D (Tgd)^ which is the correlation matrix in Eq. (15) 
associated to |T)('I'| with the restriction r G C and sgD. In [Tl] it has been shown that 
the matrices r 2 , T 3 and T 4 can be written as 


Ti = r 


AA, 


T 3 — Ti — TabiT 


-i 


B X B i 1 Bj_A, 


(29) 


and 


r 2 — m 2 y 1 m 2 , r 4 — m 2 t 3 m 2 , m 2 = 


(30) 


R, 0 
0 ~h 2 

These formulas provide an explicit representation of the matrices T, in terms of the fermion 
correlation in the finite subsystem Ai U Bi U A 2 . 

By introducing the following notation 

{..., T,... ,T',...} = Tr(... p w ... Pw’ ■ ■ ■)> (31) 


and {..., T n ,...} = Tr(... p ^ ...) as special case, from (27) we have that Tr/TJ. can be 
written as a linear combination of traces involving the matrices T*. with k G {1,2, 3,4}. 
This finally provides Trp^, which we write as follows 


T 

-L n. 


Tr Pa - YPv 


(32) 


From (21) and (27), it is straightforward to realise that T n is a combination of 4 n terms 
of the form (31) with coefficients given by integer powers of Sb 1 = {Pb^ 2 = detfTs^J. 
However, many of these 4" terms turn out to be equal when using cyclicality of the trace 
and other simple algebraic manipulations. 

In the following we write T n explicitly for 2 ^ n ^ 5, where we have computed also 
T 5 , in addition to the other ones already reported in Ref. 

• n = 2: 


t 2 = {ff} + {W r 2 } + s B , ({r|} — {r 3 , r 4 }) 


(33) 

















(35) 
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• n = 3: 

t 3 = {r?} + 3 {rf, r 2 } + 3 <j Bi ({r l5 r*} + {r 2 , r^} - 2{r ls r 4 , r 3 }); (34) 

• n — 4 : 

t 4 = {rf} + {r 1 ,r 2 ,r 1 ,r 2 } + 4 {r?,r 2 } + 2 {r?,r^} 

+ 2 s Bl ({r 1; r 3 , r 1} r 3 } + {r 4 , r 4 , r\, r 4 } + 2 {r?, r^} + 2 {rf, r*} 

+ 2 {r 4 , r 3 , r 2 , r 3 } + 4 {r 4 , r 2 , r^} - 2[ 2 {r r 3 , r 4 } 

+ {r 4 , r 3 , r l5 r 4 } + {r 4 , r 2 , r 3 , r 4 } + {r 4 , r 3 , r 2 , r 4 } + {r l5 r 2 , r 4 , r 3 }]) 

+ s 2 Bl ({r^} + 2 {r^, r^} + {r 3 , r 4 , r 3 , r 4 } - 4 {1* r 4 }) ; 

• n — 5 

n = {rf} + 5 ({rl,r 2 } + {r?,r|} + {rf,r 2 ,r„r 2 }) 

+ 5 S b , ({ry rf} + {ry ry + 2 {r 4 , ry ry + 2 {r?, r 2 , ry 
+ {r„ r 2 , r„ ry + {r„ r 2 , r„ ry + {ry r 3 , r„ ry + {ry r 3 , r 2 , r 3 } 

+ 2 {r\, r 3 , Tj, r 3 , r 2 } + {r t , r 3 , r 2 , r 3 } + {r?, r 4 , r 3 , ry + 2 {r„ r 2 , r 3 , r 2 , r 3 } 

- 2[{r=, r 2 , r 4 , r 3 } + {r?, r 2 , r 3 , r 4 } + {ry r 3 , r t , rj + {ry r 3 , r 2 , r 4 } 

+ W, r 3 , rj + {r„ r 2 , r„ r 3 , r 4 } + {r 4 , r 2 , r 4 , r 4 , r 3 } + {r 2 , r 4 , r 4 , r 4 , r 3 }]) 

+ 5 ^, ({r„ ry + {r„ ry + 2 {r 4 , ry ry + {r„ r 3 , ry ry + 2 {r„ r 3 , r 4 , r 3 , r,} 

+ {r„ r 4 , ry rj - 2[{r 4 , ry r 4 } + {r 4 , ry r 4 , r 3 } + {r 2 , r 4 , ry + {r 2 , ry r 4 , r 3 }]). 

We notice that the algebraic sum of the integer coefficients occurring in any term 
multiplying a power S Bi with p > 0 is zero. Moreover, considering only the terms which 
are not multiplied by a power 5 Bl in T n , the sum of their coefficients is 2 n ~ 1 . 


(36) 


4. Traces of integer powers of the partial transpose of the spin reduced 
density matrix 

In this section we move to the main objective of this paper which is to give a representation 
of the integer powers of the partial transpose of the spin reduced density matrix of two 
disjoint blocks with respect to kL 2 . Eisler and Zimboras in Ref. showed how to obtain 
the partial transpose of a fermionic Gaussian density matrix, a procedure which can be 


applied to the spin reduced density matrix in Eq. (19) using the linearity of the partial 

the moments of the 


transpose as we are going to show. We mention that in Ref. 
partial transpose for two adjacent intervals were studied in details using the property that 
fermionic and spin reduced density matrices are equal for this special case. 

Given a Gaussian density matrix p\y written in terms of Majorana fermions in 
A = Ai U /L 2 , the partial transposition with respect to A 2 leaves invariant the modes in 
Ai and acts only on the ones in A 2 . Furthermore, the partial transposition with respect 
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to A 2 of pa in ( fl7| ) leaves the operator P Bl unchanged (because it does not contain modes 
in A 2 ), therefore we have 


P T A = P<Tven + P B l Podd = 


T 2 


where 


To _ To i To 

P± = Peve n ± Podd> 


1 + Pl 


Bl X-2 


P+ + 


Pi 


PlpT* 


(37) 


(38) 


as clear from Eq. (19) because of the linearity of the partial transpose. The partial 
transposition of an arbitrary product of Majorana fermions A 2 (denoted shortly as 0 2 
like in the previous section) is given by the following map 


(39) 


vln) -, n r( M)o , l-/° (W mod 4) e {0,3}, 

r 2 (o 2 ) = (-1) o 2 , Tin) = | x mod 4) e {1 2}j 

where we recall that p 2 is the number of Majorana operators in 0 2 . Then, applying Eq. 

1 


(39) to (18), we find 


P T2 = 

reven 


2 ^ 1+^2 


Pol d = 


2 ^ 1+^2 


]T(-ir /2 <0 1 0 2 )0>0|, 

even 

^(_ 1 )(»-‘)/ 2 < o 1 f> Bl o 2 )o*oJ, 


(40) 


odd 


which gives the desired fermionic representation of the partial transpose of the spin 
reduced density matrix. 

At this point the moments of p T f can be obtained following the same reasoning as for 
the moments of pa • Indeed, since p± are unitarily equivalent (Pa 2 P± p a 2 — P^ 2 because 
Pa 2Peve n p a 2 = P^ven and p A 2 Podd P A 2 = ~pI, dd) and P si commutes with them, starting 
from (37) and repeating the same observations that lead to (21), one gets 

Tr(pJ)” = Tr(p£)” (41) 

Similarly to the case of the Renyi entropies considered in Sec. [3] (see Eq. ©) , the 
matrices p^ 2 are fermionic but not Gaussian. In the following we write them as sums of 
four Gaussian matrices, as done in ( [T7| ) and ( |27| ) for pa- In particular, by introducing 

1 . „+„+ _r. 1 {OiPbxO^ 


Pa = 


2 ^ 1 +fe 


E 1 “ 2 (Oi0 2 )°io\, 


Pa = 


even 

odd 


one has that the matrices in (40) become 
t 2 _ Pa + Pa 2 p\Pa 2 

reven o ’ 


2 ^ 1+^2 


Piidd ~ (PBi) 


E 

even 

odd 


Y 2 


(P Bl ) 


o f 2 ol, 


(42) 


Pa 1 ~ Pa 2 Pa 1 Pa 2 

2 i 


(43) 


telling us that p ± 2 in (38) are linear combinations of four Gaussian fermionic matrices 
occurring in the r.h.s.’s of (43). Notice that p^ en and p^ dd are Hermitian but the matrices 
defining them are not since 

(pY = PaJIPa 2 , (pf 1 )' = Pa-A'Pa,. (44) 

























Partial transpose of two disjoint blocks in XY spin chains 


11 


In order to compute the correlation matrices associated to the four matrices in Eq. 


(|43j), it is convenient to introduce 
_ / l*i 


Mo = 


0 

0 il* 2 


( 45 ) 


Then, the correlation matrices associated to p \, Pa 2 p\Pa 2 , Pa 1 and Pa 2 Pa 1 Pa 2 are given 
by 

T fe = M 2 T k M 2 , k G {1,2,3,4}. (46) 

In analogy to Eq. (32), we write the moments of pj as 


Wp?)” = M 


( 47 ) 


From Eqs. (38), (41) and (43), we have that T n is a linear combination of 4 n terms. The 
net effect is that T n can be written by taking T n and replacing F,; with Tj and 8b 1 with 
—5b 1 - The latter rule comes from the imaginary unit in the denominator of p^ Ad in Eq. 

& 

In the following we write explicitly T n for 2 ^ n ^ 5: 

• n — 2 

ft = {f ?} + {f 1, f 2 } + S B , ({f 3 , f 4 } - {f 2 

• n = 3 


3 /;) 


(48) 


T, = {p} + 3 {p, r 2 } + 34 Bl 2 ( 0 , r 4 , r 3 } - {p, p} - {r 2 , p} ; ( 49 ) 


n — 4 


= {r{} + {Ti, r 2 , r l5 r 2 } + 4 {r{, r 2 } + 2 {r{, r^} (50) 

+ 2 ( 5 / 3 , ^2{r 3 ,r 1 ,r 4 ,r 1 } + 2{r 1 ,r 2 ,r 3 ,r 4 } + 2{r 1 ,r 3 ,r 2 ,r 4 } + 2 {r 4 , r 2 , r 4 , r 3 } 

+ 4 {f 3 , f 4 , r{} - 2 {f!, f 3 , f 2, f 3} - 4 {f!, f 2 , f 2 } 

- {f l5 f 3, f 4 , f 3} - {f 4 , f 4 , f ls f 4 } - 2 {fl f ~ 2 {f l f’}) 

+ 6 2 Bl ({Ff} + 2 (r 3 , f+ {f 3, f 4 , f3, f 4 } - 4 {f 3 , r 4 }) ; 

• n — 5 

ft = {Fft + 5 ({F},F 2 } + {p,p} + {p,F 2 ,F 1 ,F 2 }) ( 51 ) 

+ 5 &Bi (2 {p,r 2 ,r 4 , r 3 } + 2{r 2 ,r 2 ,r 3 ,r 4 } + 2 {r,,r 3 ,f ,,r 4 } + 2 {r 2 , r 3 , r 2 , r 4 ) 

+ 2 {p,r 3 ,r 4 } + 2 {r 1 ,r 2 ,r 1 ,r 3 ,r 4 } + 2{r,, r 2 , r 4 ,r 1 ,r 3 } + 2 {r 1 ,r 4 ,r 1 ,r 3 ,r 2 } 

- {fi.flft- (f 1, 2 {ft, fj, f 2 } -2 {F;,F 2 ,Fi} 

- {p, p, r, : r 2 } -jp, p^p} -jp, p^p} -jp, r 3 , p, p}_ 

- 2{F 1 ,F 3 ,F 1 ,F 3 ,F 2 } - {F 4 ,F 3 ,rl,F 3 } - {rf,r 4 ,r 1 ,r 1 } - 2{p,p,f 3 ,r 2 ,r 3 }) 

4- 54 ({r„rj} + {r 2 ,r 3 } + {F„F 3 , r 2 , p} + {ft, r 4 , ft, ft} + 2 {F 1; F 3 , ft, ft, ft} 

+ 2 {ft, f l, f 4 2 } - 2 {F„ ft, ft} - 2 {ft, ft, f 4 , ft} - 2 {ft, ft, ft} - 2 {ft, F|, ft, ft}) . 
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As for T n , also in T n the algebraic sum of the integer 
multiplying a power S Bi with p > 0 vanishes. 






5. Numerical results for the ground state of the critical Ising and XX model 


The results of the previous section for the moments of the partial transpose of the reduced 
density matrix of two disjoint blocks are valid for arbitrary configurations of the XY spin 
chain: equilibrium, non-equilibrium, finite and infinite systems, critical and non-critical 
values of the parameters 7 and h. In this section we evaluate numerically these moments 
for the configurations that so far attracted most of the theoretical interest, namely the 
critical points of the XY Hamiltonian, whose scaling properties are described by conformal 
field theories. A great advantage of the present approach compared to purely numerical 
methods such as exact diagonalization or tensor networks techniques is that it allows to 
deal directly with infinite chains without any approximations, reducing the systematic 
errors in the estimates of asymptotic results. Indeed, all the numerical results presented 
in the following are obtained for infinite chains. 

We will consider two particular points of the XY Hamiltonian, namely the critical 
Ising model for 7 = h = 1 and the zero field XX spin chain (corresponding to fermions 
at half-filling) obtained for 7 = h = 0. The scaling limit of the former is the Ising CFT 
with central charge c = 1 / 2 , while the scaling limit of the latter is a compactified boson 
at the Dirac point with c = 1. 

The CFT predictions for the moments of both reduced density matrix and its partial 
transpose have been derived in a series of manuscripts and they are reviewed in Appendix 
[A] For both models we consider the case of two disjoint blocks of equal length £ embedded 
in an infinite chain and placed at distance r. We numerically evaluate the moments of pA 
and p T f using the trace formulas of the previous sections for n = 2,3,4, 5 and we compute 


the ratio R n (defined in Eq. (11)): 


_tv( p JT 

■tin = 


Tip 


(52) 


whose (unknown) analytic continuation for n e —» 1 would give the negativity. Notice that 
from Eqs. (32) and (|4T[) we have that R n = T n /T n . In the scaling limit (i.e. £ 1 r —>• 00 


with ratio fixed) the ratio R n converges to the CFT prediction (cf. Eq. (81) in Appendix) 
written in terms of the four-point ratio x, which is 

2 

(53) 


X = 


+ r 


when specialised to the case of two intervals of equal length £ at distance r. 


5.1. The critical Ising chain 

The negativity and the moments of p 1 / for the critical Ising chain in a transverse field 
have been already numerically considered in Ref. [T5] by using a tree tensor network 
algorithm and in Ref. [19] by Monte Carlo simulations of the two-dimensional classical 
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Figure 2. The ratio R n between the integer moments of pa and pfj for two disjoint 
blocks of length £ at distance r embedded in an infinite critical Ising chain. We report 
the results for n = 3,4, 5 as function of the four-point ratio x for various values of £ (and 
correspondingly of r). For large £, the data approach the CFT predictions (solid lines). 
The extrapolations to £ —> oo -done using the scaling form (541 are shown as crosses 
and they perfectly agree with the CFT curves for n = 3 and 4, while for n = 5 the fits 
are unstable and the extrapolations are not shown. The last panel shows explicitly the 
extrapolating functions for two values of x and n = 3,4. 


problem in the same universality class. However, the finiteness of the chain length did not 
allow to obtain very precise extrapolations to the scaling theory for all values of n and of 
the four-point ratio x. We found, as generally proved [9], that i?2 is identically equal to 1. 
In Fig. [2] we report the obtained values of R n for n — 3,4, 5 as function of x for different 
values of l. It is evident that increasing £ the data approach the CFT predictions (the 
solid curves). We can also perform an accurate scaling analysis to show that indeed the 
data converge to the CFT results when the corrections to the scaling are properly taken 
into account. 

It has been argued on the basis of the general CFT arguments [41], and shown 
explicitly in few examples [35] 05] 06] both analytically and numerically, that Tr p n A displays 
‘unusual’ corrections to the scaling which, at the leading order, are governed by the 
unusual exponent S n = 2 h/n where h is the smallest scaling dimension of a relevant 
operator which is inserted locally at the branch point [53]. For the Ising model it has 
been found that, in the case of two intervals, h = 1/2 [ 351 PfT j. From the general CFT 
arguments in Ref. m, we expect the same corrections to be present for Tr(pj) n because 
they are only due to the conical singularities. Unfortunately, the corrections to the scaling 
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in Fig. [2] cannot be captured by a single term, because subleading corrections become 
more and more important when n increases, as already pointed out in Ref. |T8]. Indeed, 
corrections of the form £~ m / n for any integer m are know to be present [351 HU[57]. Thus 
the most general finite-^ ansatz is of the form 


JD D 1 

ll n Tt n 


CFT 


Or) + 


Tn ] I 


£R 


x) r0 2) 


IP' 


x) r j 3) (x) 


£3/r, 


(54) 


The variables r n (x) are used as fitting parameters in the extrapolation procedure. The 
number of terms that we should keep in order to have a stable fit depends both on n and on 
x. For each case we keep a number of terms such that the extrapolated value at £ —* oo is 
stable. In any case we never keep corrections beyond the order 0(f' _1 ). The results of this 
extrapolation procedure for n = 3 and 4 are explicitly reported in Fig. [2] The agreement 
of the extrapolations with the CFT predictions is really excellent, at an unprecedented 
precision compared with fully numerical computations [TH, THj- Conversely, we find that 
for n = 5 the extrapolations are still unstable because of the large number of terms we 
should keep in order to have a precise enough extrapolation. 


5.2. The XX chain 


We now move to the study of the powers of frf for the XX model in zero held. There are 
no previous numerical studies of this paradigmatic model. We again consider the ratios 
R n for n = 2,3,4, 5 and we again find that R 2 is identically equal to 1, as it should be. 
In Fig. [3] we report the obtained values of R n for n — 3,4, 5 as function of x for different 
values of £. It is evident that increasing l the data approach the CFT predictions (the solid 
curves). We should however mention a very remarkable property. It has been observed 
that Tr p n A shows oscillating corrections to the scaling |35L IT71141 j . which for zero magnetic 
held, are of the form (—l)h These oscillations however cancel in the ratio R n and the 
corrections to the scaling are monotonous, a property which makes the extrapolation to 
infinite £ slightly simpler. 

Also in this case we can perform an accurate scaling analysis to show how the data 
converge to the CFT results when the corrections to the scaling are properly taken into 
account. For the XX model, the leading correction to the scaling is governed by an 
exponent 5 n = 2/n, which means that they are less severe than in the case of the Ising 
model as it is also qualitatively clear from the figure. We then use the general finite- 
ansatz 


JD _ D 1 

1 t n 1 l n 


CFT, 


x) + 


Tn\x) 

£2/n 


+ 


r ( n\ 


£R' 


x) r j 3) (x) 


£6/r, 


(55) 


and, as in the case of the Ising model, we keep a number of fitting parameters which 
make stable the extrapolation at £ —> 00 . The results of this procedure for n = 3 and 
4 are explicitly reported in Fig. [3] The agreement of the extrapolations with the CFT 
predictions is excellent. Also for the XX chain we find that for n — 5 the fits are unstable. 
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Figure 3. The ratio R n between the integer moments of pa and pff for two disjoint 
blocks of length l at distance r embedded in an infinite XX chain at zero field. We report 
the results for n = 3,4, 5 as function of the four-point ratio x for various values of £ (and 
correspondingly of r). For large £, the data approach the CFT predictions (solid lines). 
The extrapolations to £ —> oo -done using the scaling form (551 are shown as crosses 
which perfectly agree with the CFT curves for n = 3 and 4, while for n = 5 the fits are 
unstable. The last panel shows explicitly the extrapolating functions for two values of x 
and n = 3,4. 


6. Two disjoint intervals for free fermions 


In this section we consider the partial transposition for two disjoint blocks in the fermionic 
variables. This problem was already addressed by Eisler and Zimboras (3?J, but a detailed 
numerical analysis was not presented. For fermionic variables there is no string in B\ 


connecting the two blocks (cf. Eq. (16)). Thus the partial transpose of fermions can 
be obtained from the formulas derived in the previous sections by discarding the string 


of Majorana operators (16), i.e. by replacing Pb 1 with 1. Performing this replacement, 


many simplifications occur in the formulas found in Sec. [3] and Sec. [4] as we will discuss 
in the following. 


6.1. Renyi entropies 

By definition the fermionic reduced density matrix is Gaussian with correlation matrix 
Ti defined in the Sec. [3j i.e. 

Ti Va = {U}- 


(56) 
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It is however instructive to recover this result from the formulas in Sec. |3] in order to set 
up the calculation for the partial transpose. 

Making the replacement Pb 1 —> 1, the reduced density matrix of the two disjoint 
blocks given in Eqs. © and ( [ 18 ] ) becomes 

Pa = Peven + Podd; (57) 

whereJU 


1 


Peven 


2^i+4 Z —i 


( 0 , 02)0101 


F _ 
Podd 


2 ^ 1+^2 

odd 


( 0 , 0 2 ) 0 \o\. 


(58) 


Moreover, p ^ 1 defined in Eq. ( 23 ) is replaced as p B jf —>• pa and therefore, from Eqs. ( 27 ) 
and ( 57 ) we conclude that p B = p\. As for the correlation matrices r ?; , since p Bl —>■ p A , 


2 - 


it is obvious that T 3 —> Ti and T 4 —> T 

Summarising, we conclude that the fermionic Tr p n A is found by making in Eq. (32) 
the following replacements 

<5bi —> 1, r 3 y ri, T4 —)■ T 2 - (59) 

Performing these substitutions in the explicit examples given in Sec. [3] for 2 ^ n ^ 5, 


it is straightforward to find Eq. (56), which is just the obvious result that the fermionic 


density matrix is the Gaussian operator with correlation matrix given by Ti. 

As a further check of our numerical codes, we numerically calculated Trp^ using Eq. 


(56) (as was already done in Ref. H), obtaining that on the critical lines in the scaling 
limit it converges to 

Tr p n A ra „ xl2A , (60) 


[£ 1 ^ 2(1 - £)] 


that corresponds to J- n (x) = 1 identically in the general CFT formula (73). Indeed this 
result was already proven in the continuum free fermion theory m ■ 


6.2. Traces of integer powers of the partial transpose 

We are now ready to set up the formulas for the moments of the partial transpose, as 
already derived by Eisler and Zimboras |37] , but numerically studied only for the case of 
adjacent intervals. 

Once again, Tr(pj) n in the fermionic variables is obtained by replacing Pb, with 
1 in the formulas reported in Sec. [4j Performing this replacement in Eq. (37) we get 


Pa = Peven + (Po< dd ) r2 = P+ an d in Eq. (42) it gives p ^ 1 — » p\. These observations 


together with Eq. M3|) lead to 


P T 1 


1 -i 


Pa o Pa-zPaPaii 


(61) 


2 2 

which is exactly the same result obtained in [37] (for direct comparison, set (0 + ) t here = 

(/bOhere Slid (0_)there {PA 2 PAPA 2 ) here)• 


§ To compare our notation with the one used in m , Set (p+) there — (Peven)here &nd (p_) there — (Po F dd) here* 
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Figure 4. The ratio R n between the integer moments of pa and pff for two disjoint 
intervals of length £ at distance r for the tight-binding model at half-filling. We report 
the results for n = 3,4,5 as function of the four-point ratio x for various values of £ 
(and correspondingly of r). For large £, the data approach the CFT predictions (solid 
lines). The extrapolations to £ —> oo -done using the scaling form (55)- are shown as 
crosses and they perfectly agree with the CFT curves for n = 3,4,5. The last panel 
shows explicitly the extrapolating functions for one value of x and n = 3,4, 5. 


The last remaining step is just to write these formulas in terms of correlation matrices. 
Given that p^f —> p \, it follows that we should perform the replacements T 3 —> IV and 
T 4 —^ r 2 in order to get the moments of the partial transpose in terms of the correlation 
matrices. Summarising, the fermionic Tr(p^ 2 ) n are given by the formulas in Sec. [ 4 ] 
performing the replacements 


5b 1 —>• 1, T 3 Ti, T4 —> IV 


(62) 


Writing Tr(/y 4 2 ) n = Vf /2 n , and performing the replacements in the formulas for 

2 ^ n ^ 5 given in Eqs. pjj ), (49), (50) and Q, we find 


Tf = 2{r 1 ,r 2 }, ~ _ (63) 

tJ= -2{r?} + 6{r?,r 2 } 1 ~ _ (64) 

tJ = -4{r^} + 4 {ri,r2,r 1 ,r 2 } + 8{r?,r2}, _ (65) 

T( = - 4 {rf} - 20 {Ti r 2 } + 20 {rf, r^} + 20 {r ?, r 2 , r 1? r 2 }. (66) 


Notice that the final expressions are very compact compared to the much more 
cumbersome spin counterparts. 
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We are now going to evaluate numerically the moments of the reduced density matrix and 
its partial transpose. We can study the problem for arbitrary values of h and 7 entering 
in Eq. ([3]), but in the following we focus on the most physically relevant fermionic system 
with h — 7 = 0 , i.e. the tight binding model 


1 L 

H=~Y 

2 ^ 

i —1 


CjCi+l + 4+l C i 


(67) 


at half filling (kp = 7 r/ 2 ). In the scaling limit, the tight binding model is described by a 
CFT with c = 1. 


The numerical results for the ratios R n = Tr(p^ 2 ) n /(Trp^) are reported in Fig. [ 4 ] 


as 


function of the four-point ratio x for different £ and for n = 3,4,5 (we checked that R 2 = 1 
identically, as it should). We also derived asymptotic CFT predictions for the fermionic 
moments of the partial transpose, but their derivation is too cumbersome and beyond the 
goals of this manuscript. We will report the derivation in a forthcoming publication [48 J 
and we limit here to give the final results for n = 2, 3,4, 5. In order to have manageable 
formulas we introduce the shorts 


2e 

L 2S 


|0[e](O|r(x))r 


( 68 ) 


|@(0|r(x))r 

where 0 is the Riemann Theta function defined in Appendix [A} In terms of the 0 function 
the fermionic ratios R n are given by 


2 Ro 


(1 - x) 4A2 

4 R 3 

(1 - X) 4A 3 

8R4 

(1 — x) iAi 
16 R 5 

(1 - :r) 4A5 


0 

L 1 J 


-2 + 6 


-4 + 


4 + 20 


0 

0 

0 

0 

r 0 
0 


0 

0 


+ 4 


0 1 
1 

+ 20 


0 

1 

0 

0 


20 


0 

0 


(69) 


(70) 


(71) 


(72) 


where the matrix f has been defined in Eq. (80) and the exponent A n in App. A 

It is evident from Fig. [4] that the lattice numerical results approach the CFT 
predictions depicted as solid lines for all n. As in the spin case, we can perform a careful 
finite i analysis to take into account corrections to the scaling. The leading correction is 
expected to be of the form £~ 2 / n and subleading ones to be integer powers of the leading 
one. The finite £ ansatz is then given by Eq. (55) and again to have an accurate description 
of the data we keep a number of fitting parameters which make stable the extrapolation at 
£ —> 00. The results of this extrapolation procedure for n = 3,4, 5 are explicitly reported 
in Fig. [4} The agreement of the extrapolations with the CFT predictions is excellent also 
for n = 5 as a difference compared to the spin counterpart. 
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7. Conclusions 

We have shown that the partial transpose of the reduced density matrix of two disjoint 
spin blocks in the XY spin chain can be written as a linear combination of four 
Gaussian fermionic operators, fully specified by their correlation matrices (denoted as Ik, 
i — 1, 2, 3,4 in the text) which have been explicitly calculated in terms of the correlation 
matrix of the subsystem formed by the two blocks joined with the finite part between 
them. This construction allows to calculate the moments of the partial transpose in 
generic configurations of the spin chain. In this manuscript we focused on the ground 
state of Ising and XY chain, but the approach is more general and can be used for 
arbitrary excited states, thermal density matrices, non-equilibrium situations etc. 

The obtained representations of the moments of the partial transpose allow us to 
study in an exact manner infinite chains and very large subsystems, drastically reducing 
the systematic errors in the approach to the scaling limit. We found that for the ground 
state of the critical models the moments of the partial transpose agree (with high accuracy) 
with the recent CFT predictions after the corrections to the scaling are properly taken 
into account. We also studied numerically the integer powers of the partial transpose in 
the fermionic degrees of freedom (described in Ref. [37], but not numerically studied). 
Even in this case we find that the moments agree perfectly with the CFT predictions in 

US- 

The main open problem left by this manuscript for two disjoint blocks of a spin 
chain is whether it is somehow possible to obtain the negativity from the correlation 
matrix (the problem is also present for the fermionic degrees of freedom [37]). A similar 
problem is also open for the entanglement entropy since integer moments are obtained 
in a similar fashion [41], but one has no access to the spectrum of the reduced density 
matrix and hence to the entanglement entropy. From the practical point of view, it 
has been recently shown that if one knows a relative large number of integer moments, 
rational interpolations provide accurate estimates of the analytic continuations [ 491130] 
(and hence of entanglement entropy and negativity). However, a deeper understanding 
of these analytic continuations would be highly desirable. 
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A. CFT results for entanglement entropy and negativity of two disjoint 
intervals 


The moments of the reduced density matrix of two disjoint intervals for CFTs have been 
studied in a series of manuscripts [HI EH [52], [53J Ell EH EE EH ESI EH1 EOj EH E3 [63] . 
These results have been derived using earlier findings for the partition functions of CFTs 
on Riemann surfaces with non vanishing genus |64| . In this appendix we review the 
main results (especially from Refs. pH EH EEj) which are useful for the comparison 
with numerical results. We mention that some universal results are also known in higher 
dimensions both from field theory [65] and holography [63, 661. f671j. 

From global conformal invariance, we know that Trp^ for two disjoint intervals admits 
the general scaling form (choosing, without loss of generality, the endpoints of the intervals 
in the order u\ < Vi < U 2 < v 2 )\ 


Tr p n A 


( (u 2 - Ul ){v 2 - v x ) \ " 

\(V 1 -U 1 )(V2-U 2 ){V2-U 1 )(U2-V 1 )) 


(73) 


where A n 
ratio 


c{n — l/n)/12, being c the central charge. The variable x is the four-point 


{ui - Vi){u 2 - v 2 ) 

(■ ui - u 2 )(v 1 - v 2 y 


(74) 


Given the order of the points we have 0 < x < 1. The prefactor c n is non-universal, 
but can be exactly fixed from the exact calculation of the entanglement entropy of one 
interval. 

The difficult task of CFT is to have an exact representation for the universal function 
X n {x) normalised so that 7-^(0) = 1. This universal function has been analytically derived 
for the compactified boson (with central change c = 1) [521153] and for the Ising CFT (with 
c = 1/2) [561157] . as well as for other conformal theories which however are not of interest 
for this paper. Concerning the compactified boson, we are only interested in the value of 
the compactfication ratio corresponding to the scaling limit of the XX spin chain which 
is the so called Dirac point. For the Ising CFT and at the Dirac point (which describe 
respectively the scaling limit of the critical Ising chain and the critical XX model), the 
function lF n (x) reads [55] 


Rising / 7 ,N = Ee MgKlMH 7rDirac / \ = Ee I'® [ e ] ( t ( x ) ) |~’ 

n M 2 n_1 |0 (t(x))[ n 2 n ~ 1 |0(r(a;))| 2 

where 0[e](f2) is the Riemann theta function, which is defined as follows [68] 


(75) 


i7r(m+e) t -r2-(m+e)+27ri 

m G Z n— 1 

being Ha(n»l)x(n- 1) symmetric complex matrix with positive imaginary part and 
e is the characteristic of the Riemann theta function, which is defined by a pair of n — 1 


@[ e K^) = e 


e = 


‘ £ ' 


- £ 1 ,. 

• 5 ^n—1 

. s. 


L <5i,. 

• ? On—1 - 


. (76) 
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dimensional vectors made by e i: 8i G {0,1/2}. In Eq. (76) we have to sum over all the 
characteristics e. The elements of the matrix t(x) in Eq. (75) read 


n— 1 


T(x) rs = i — > sin(7rfcm) 
n 


fc=i 


■ 2 Fi(k/n, 1 — kjn\ 1; 1 — x) 
2 F 1 (k/n, 1 - k/n ; 1; x) 


cos 


k . 
27T —(r — s) 
n 


(77) 


where x € (0,1) and 2 F\ is the hypergeometric function. 

Also the moments of the partial transpose correspond to a four-point function of 
twist fields in which two of them have been interchanged [9]. Consequently also these 
moments admit the universal scaling form 


Tr(p' 


T 2 \n 
A ) 


= c: 


(u 2 - Ui)(v 2 - Ui) 


(ni - Ui)(v 2 - U 2 )(v 2 - Ui)(u 2 - Vi) 


2A„ 


Qn(x) 


(78) 


with c n the same non-universal constant appearing in Eq. (73) and G n {x) a new universal 
scaling function. Exploiting the fact that the above moments correspond to the exchange 
of two twist fields, it has been shown that G n (x) and F n {x) are related as pun! 


Qn{x) = (1 - x) 4A " J7 


X 


X — 1 


(79) 


but some care is needed to take the analytic continuation of the function F n (y) to negative 
argument y (see for details to. This result is equivalent to say that Tr(pj) n is given by 
Eq. (75) in which the period matrix t(x) is replaced by 


t(x) —> rlx) = T 


x 


x — 1 


Thus, the CFT prediction for ratio in Eq. GD 

4A„ T„(X/(X - 1)) 


IS 


R 


CFT 


(d) = (1 — x) 


FJx) 


( 80 ) 


(81) 


in which the universal constants c n as well as the dimensional part of the traces canceled 
out leaving a universal scale invariant quantity. Thus, in order to study this quantity 
we do not need an a priori knowledge of the constants c n which anyhow are known both 
for the XX [32] and the Ising [33( [69] spin chains. It is worth recalling that the analytic 


continuation to non-integer n of the ratios (81) for two intervals are not yet known even for 
the simpler cases and consequently also the negativity is eluding an analytic description. 
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